Far-field image of Veselago lens 

C. Y. Li, J. M. Holt, and A. L. EfroJ] 

University of Utah, Salt Lake City UT, 84112 USA 
(Dated: February 2, 2008) 

Abstract 

It is shown that perfect imaging of a point source both in near- and far-field regions contradicts 
electrodynamics although "superlensing" is impossible only in the far-field region. These general 
statements are illustrated by detailed study of evanescent wave propagation in a photonic crystal 
that is known to be a left-handed medium. An analytical expression for the intensity distribution 
near the far-field focus is obtained. This distribution contains some novel features and establishes 
a new "diffraction limit" for flat lenses. It is generalized for multiple sources located at different 
points. The distribution is in very good agreement with computer simulation. 

PACS numbers: 78.20. Ci,41.20.Jb, 42.25.-p 
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I. INTRODUCTION 



The left-handed medium (LHM), defined by Veselago^ as a medium with simultaneously 
negative and real n and e, recently has attracted much attention mostly because of the 
negative refraction at its interface with a regular medium (RM). This effect allows creation 
of a unique device called the "Veselago lens" . This lens is a slab of LHM inside a RM with a 
condition that both media have the same isotropic refractive index and the same impedance. 
Interest in LHM's significantly increased after the work by Pendry^ which argued that the 
Veselago lens is a "perfect lens" in the sense that it gives a perfect image of the point source. 
This statement is based upon the observation that the evanescent waves (EW's) of a form 
exp(ik y y — kx) that usually decay in the near-field region are amplified by the LHM. Pendry 
claimed that the amplified EW's restore a perfect image in both the near-field and far-field 
regions. Following Veselago, Pendry considered a hypothetical LHM (HLHM) with negative 
and real fi and e. He did not present a solution in coordinate space. This was done by 
Ziolkowski and Heyman^ and their solution reveals a fundamental problem with Pendry's 
idea. The solution diverges exponentially at each point of a 3-D domain near the focus, 
just where the fields of the EW's increase due to amplification. This was found out almost 
simultaneously by three groups4*^ Pokrovsky and Efros 5 used diffraction theory, known to 
be exact at small deviations from geometrical optics, and found a finite width of the focus 
in the far-field region. Haldane^ argued that the problem in general does not have a solution 
in the framework of macroscopic electrodynamics. 

In this paper, we study the far-field imaging of a two dimensional Veselago lens. In the 
first section, we present two-dimemsional diffraction theory and use it to find an analytical 
expression for the intensity near the focal point of the Veselago lens in the HLHM model. 
As in the three-dimensional case, it is exact in the far-field region and does not suffer from 
any of the problems introduced by EW's. We present also a general proof that the perfect 
imaging of a point source contradicts to wave optics. All results obtained in this section are 
independent of the microscopic structure of the LHM that forms the lens. 

In the second section, we consider a uniaxial photonic crystal (PC) and present computer 
simulations to show that, in the frequency range near the T-point of the second Brillouin 
zone, it behaves as a real LHM, described by Veselago, with respect to propagating mode^ 
but it does not provide universal amplification of EW's. The reason is that, for the mathe- 
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matical description of the EW in the PC, one should take into account the spatial dispersion 
(k-dependence) of p and e. However, the amplification may (or may not) appear due to some 
other reasons rather than negative p and e. For example, the amplification might be due 
to surface waves that are not related to the dielectric properties of the PC. Such an ampli- 
fication may provide an improvement of the image in the near-field region, but it does not 
affect the image near the far-field focal point of the Veselago lens. Using these results, we 
provide a wave theory, more general than diffraction theory, which gives exactly the same 
analytical expression for the field intensity near the far-field focus as diffraction theory. 

In the third section, we present the results of a computer simulation for focusing of both a 
thin (near-field regime) and thick (far-field regime) Veselago lens and show that the far-field 
image is described by the analytical theory without the contribution of EWs with very good 
accuracy, while the near-field image may have "superlensing" features. We also define a class 
of field distributions which gives a perfect far-field image without any "superlensing". In 
fact, the Fourier transform of such fields should not include the EWs at a given frequency. 

II. TWO-DIMENSIONAL DIFFRACTION THEORY AND IMAGE OF THE 
VESELAGO LENS 

The two dimensional Green function of the Helmholtz equation in a RM can be written 
in a form Gr = (i / 4) H^\pk ) , where k = ujn/c, the Hankel function = J +iN , and 
p = y/Jx — x ) 2 + (y — Do) 2 , where x and y are the coordinates of the point source. The 
function G R oscillates as exp ik p at large p. The Green function in the LHM should have a 
form Gl = — {i/ 4)Hlp* (pk®). At large p, it oscillates as exp—ik p, while at small p, it has 
the same form as Gr. 

We consider a two dimensional problem in the plane x — y with a point source at a 
distance a from the LHM slab (see Fig. 1). It creates a magnetic field in the z-direction. 
We want to find the field near both foci. It is easy to show that the perfect image of this 
source is impossible. Indeed, if the image has the same form as the source, the magnetic field 
near the focus would obey the inhomogeneous Helmholtz equation with a delta-function at 
the focus, while the focal point would not contain any source. Thus, this solution does not 
obey the homogeneous Helmholtz equation near the focal point. Note that this argument is 
general and does not specify the nature of the lens. 
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FIG. 1: Refraction of light outgoing from a point source at x = —a and passing through a slab of 
LHM at < x < d described by anomalous Snell's law. A second focus (the image) is formed at 
x = 2d — a. The arrows show the direction of the wave vector. 

To find the distribution of the field, we use diffraction theory based upon Huygen's 
principle.- The theory is exact at small wavelength when deviation from geometrical optics 
is small. Suppose that the field H(y) is given on a line x = 0. Since the wavelength is 
small, the field can be represented locally as a plane wave with wave vector (k x , k y ). The 
field inside the LHM can be found by implementation of Huygens's principle in the two- 
dimensional case: 



Here 8 is the angle between the direction of the ray and the x-axis. 

To find b r i, we use the same method as Ref. 9. Assume that the field to the left of the 
LHM slab has a form of a plane wave, namely, H = Hoexpikox. Then Eq. (1) should give 
us H = H exp —ikox inside the LHM slab if k x 3> 1. From this condition one can find b r i. 
Since we assume k$x 3> 1, the Hankel function can be taken in the form 



One gets cos(#) = I for a plane wave at normal incidence. To perform the integration in Eq. 
(I) we can expand the exponent with respect to y' — y and neglect y' — y in denominator. 
The result is b r i = 2ki. 




(1) 




(2) 
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Now we find the field in the LHM created by the point source at x — —a, y = (see 
Fig. 1). The field of the source at x < has a form H s (x,y) = iH HQ 1 \k p), where H is 
the amplitude and p = yj{x~+ a) 2 + y 2 . Thus we should substitute H(y') = H s (0,y') and 



cos(6) = a/ \J a? + y' 2 into Eq. (1). Then inside the LHM slab (0 < x < d) the field has a 
form 

/(l poo exp ik (Va 2 + y' 2 - y 1 x 2 + (y - y') 2 ) 



H(x,y) = H — ^ TJ r—^dy'. (3) 

W-oo (a 2 + y' 2 )^ [x 2 + (y- y') 2 } 1/4 

If one considers the field near the first focal point (y = 0, x = a) by assuming y <^ a, 
\x — a\ <^ a, ka ^> 1 and introducing rji = ko(x — a) and ( = koy one obtains 

H(x, y) = H - j ^ ^_ _ dt. (4) 

To find the field in the RM region x > d we can use the Huygens's principle in the form 

/oo 
H(y')G R (k ^x 2 + (y'-y) 2 ) cos(9)dy', (5) 
-oo 

where integration is performed over the right interface. Using a plane wave, as done above, 
one finds that the coefficients b r i = —h\ r . The field H(y') at the right interface can be 
obtained from Eq. (3), by setting x = d. In the vicinity of the second focus (x = 2d — a and 
y — 0) one gets 



tt( \ tt f expt(Ct + yVl-t 2 ) 

H(x,y) = H - dt, (6) 

where r\ = k (x — 2d + a). 

The field in the y-direction at x = 2d — a is 

H{2d -a,()= H - f ^2=dt = iH J {Q (7) 

and the horizontal profile of the second focus at y = is 

2i f 1 \cos(riq) + isin(riq)] , ^ . . T . 

77, 0) = if„- / 1 = — rfg = #o -Ho fa + ^o(77)] , (8) 

^ Jo V 1 - g 2 

where Jo (z) and H (-2) are the Bessel and the Struve functions, respectively. The dimen- 
sionless intensity \H(x,y)/ Hq\ 2 near the foci is shown in Fig. 2 in both perpendicular (x) 
and lateral (y) directions as a function of a distance in units of wavelength A = 2-w/ko. 

The intensity near the focus is not symmetric with respect to x and y: the width of the 
distribution in the direction perpendicular to the slab is larger than in the lateral direction 
though the source field depends only on p. 




FIG. 2: Distribution of dimensionless magnetic energy near the foci of the Veselago lens as a 
function of £ and r/ as given by Eq. 6. 

Another important feature of this imaging is that the image in the lateral direction is 
H(2d — a,y) = iH J Q (\y\k ) exp — iut. If we consider real fields, the part of the source field 
that is proportional to sm(ut), namely HoJo(pko) sin(a;t), has a perfect image in the lateral 
direction while the part of the source with cos (cut) has a zero image in this direction. This 
follows only from the flat (non-axial) geometry of the lens that is assumed to be infinite in the 
lateral direction. The perfect lateral imaging of the Bessel function Jo does not contradict to 
general principles because this Bessel function obeys the homogeneous Helmholtz equation. 

III. PROPAGATION OF EVANESCENT WAVES THROUGH A LEFT-HANDED 
MATERIAL MADE OF PHOTONIC CRYSTAL 

We consider a dielectric photonic crystal (PC) that is known to be a LHM in some 
frequency range. 7,8 Here we show that its electrodynamics differs from that of the HLHM, 
namely that the negative [i and e of propagating waves do not provide amplification of EW's. 

A large group of works pioneered by Notomiifi and Luo et al . 11 ^ 12 claims observation of 
negative refraction, focusing, and even superlensing in a PC due to photons with wave vectors 
deep in the Brillouin zone. Propagation of such photons can be described macroscopically 
only by taking into account spatial dispersion (k-dependence of e and //). The theory of 
Veselago can hardly be generalized for the case of spatial dispersion mostly because of an 
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FIG. 3: Five lowest bands of the photonic spectrum of a 2-D PC with a period d. The unit cell of 
the PC is shown in the inset. The PC consists of a square lattice of circular cylindrical air holes 
in a dielectric matrix with e m = 12, fi m = 1. The radii of the holes R = 0.35d. 

extra term in the definition of Poynting's vector.^ Following Ref. 11, we do not consider 
such a medium left-handed and do not discuss it here. 

It was shown recently,^ that a 2-D dielectric uniaxial PC made of non-magnetic materials 
can behave as a LHM with negative e and \x if it has a negative group velocity in the vicinity 
of the r-point. This was proved for propagating modes only. Experimental demonstration 
of negative refraction in a metallic PC using the modes near the r-point has been presented 
in Ref. 14. 

In this paper, we consider p-polarization for both EW's and propagating modes in a 
uniaxial PC. Fig. 3 shows the elementary cell in the plane x — y and the spectrum of p- 
polarized propagating waves. The dashed line shows our working frequency (u)d/27rc = 0.33) 
for all data below. Computations are performed using the method of finite elements with 
the software "FEMLAB". 

From the same numerical solution of the microscopic Maxwell's equation that gives the 
spectrum, one gets microscopic fields h z (x, y), e y (x, y), and e x (x, y) corresponding to a given 
frequency of propagating modes. The microscopic fields are the Bloch functions while the 
macroscopic fields have a form B z =<h z >, E =<e>, where <...> means averaging over the 
unit cell. Near the r-point, the Bloch functions have a small k and the macroscopic fields 
of the propagating modes are plane waves that obey the macroscopic Maxwell equations. 
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As follows from Ref.—, the values of \i zz = fi and e xx = e yy = e that describe propagation 
of these waves are negative because group velocity is negative. Here, z is the axis of the 
crystal. It has been shown previously^ that the values of e and fi can be expressed through 
the integrals of Bloch's function at k — 0. It is important to understand that, in a system 
with spatial dispersion, magnetization and displacement currents are inseparable- 13 Our 
parameters e(u) and fi(uj) are not exclusive properties of the medium at a given frequency, 
as in macroscopic electrodynamics without spatial dispersion. Rather, they are properties 
of a particular mode of this medium, whose fields are used for the calculations. We argue 
below that they are not applicable to the EW's. The theorem^ connecting the signs of both 
e(u) and fi{ui) with the sign of group velocity is also not applicable to the EW's. Near the 
T-point n ~ fc 2 .— Since the dispersion law u(k) is also isotropic, one obtains fi that is a 
function of uj only. This approach is valid for an infinite PC, but may create problems with 
boundary conditions— because, in coordinate-space, k 2 should be considered as a Laplacian 
operator. We will refer to this problem later on. For EW's, this approach fails completely 
because their dispersion law is anisotropic. 

To demonstrate the Veselago lens, one should know fi and e of the RM medium around 
the PC. In this paper, we find the value of the refractive index n = ^Jejl from the function 
oo(k), shown in Fig. 3, using the definition u = kc/n. The second condition that provides 
matching of impedance is the absence of reflection of the incident plane wave in a wide range 
of incident angles. 

Using these two conditions, we get that, at the working frequency, the RM material 
should have 

e = 1.125, fi = 0.08. (9) 

The physical reasons for the appearance of fi in a non-magnetic PC are discussed in Ref. 
7. If not stated otherwise, we use these values for the homogeneous regions surrounding the 
PC slab in the computations below. 

Fig. 4 shows the result of our simulation of light propagation through a PC slab sur- 
rounded by a homogeneous medium. One can see the negative refraction of light coming 
in and going out of the PC with equal absolute values of incident and refracted angles and 
without any visible reflections from the interfaces. We verified this condition over a wide 
range of the incident angles (from to 60 degrees). An animation^ shows that the wave 
fronts are moving to the right outside the PC slab and to the left inside the slab. All of 

8 



FIG. 4: Negative refraction of light propagating through a PC slab surrounded by a homogeneous 
medium at the working frequency. Periodic boundary conditions are used in the vertical direction. 
The arrows show the direction of the wave vector. 

these phenomena exactly correspond to Veselago's theory.— 

Now we consider EW's in the same PC slab embedded in the same RM as in Fig. 4. Our 
results are shown in Fig. 5. The boundary conditions implied for an incident wave at x = 
are h z = exp(ik y y) and dh z /dx = —nh z . Note that, due to reflection, the total field shown 
in Fig. 5 may not coincide at x = with the incident one. To obey Maxwell's equation, 
we imply the dispersion law of the EW's in a form k 2 — k 2 = u 2 n 2 /c 2 . The distance from 
the left boundary of the RM (point x = 0) to the PC slab is 1/2 of the slab thickness. In 
Fig. 5(a) the PC slab is substituted by the HLHM that has negative e and \x with the same 
absolute values as in the RM. One can see a strong amplification of the EW first discovered 
by Pendry. 2 However, the PC slab gives a very different picture. First, there are oscillations 
inside the PC with the period d of the PC. The result for the envelope strongly depends on 
the way the surface is cut. Our surfaces are perpendicular to the [10] direction, but they are 
cut either across the holes (AH) or between the holes (BH) (see insert in Fig. 5 (b),(c)). 

In the case of the BH surface (dashed lines), the EW decays in the PC slab in the same way 
as in the RM without any signature of amplification. However, there is amplification at the 
left AH surface for the thick slab and at both AH surfaces for the thin slab. We believe that 
these amplifications are due to surface waves (SW's), excited by EW's. Since the picture is 
not symmetric with respect to left and right interfaces, this excitation is not resonant. These 
SW's have nothing in common with left-handed properties or with the surface polaritons in 
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FIG. 5: Magnetic field of the EW in a HLHM (a) and in a PC slab (b),(c) surrounded by a 
regular homogeneous medium at the following parameters for the EW: k y = (\/ r 5/2)ko, k = 0.5ko, 
ko = ujn/c, n = ^feji = 0.3. The cross section y = that goes across the holes is shown. The 
results are similar if the cross-section y = goes between the holes. The thickness of the PC slab 
is 20 d ~ 2 A (b) and lOd ~ A (c). Here A = 2^ jk^ is the wavelength. The solid lines are for the AH 
surfaces, the dashed lines are for the BH surfaces. The period of oscillation is d. 



the HLHM, considered by Ruppin^ because they depend on the properties of the surface, 
while the polaritons should have a wavelength larger than the lattice constant, so that their 
spectra depend on the e and \i of the bulk materials. We think these SW's appear solely 
due to the termination of the periodic structure. They have been studied before by many 
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FIG. 6: An EW propagating in a direction not parallel to the surface of the PC slab exhibits 
negative refraction at the interface. The arrows show the directions of the propagating wave 
vector. The EW decays in the plane x — y in the direction perpendicular to the propagating wave 
vector. 

author o 12 i 17 i 18 with and without connection to focusing. 

Thus, we have found that EW's cannot be described in terms of the same negative e and 
/j, that are found for the propagating waves. The formal reason is that the dispersion law for 
the EW is anisotropic even near the T-point. Therefore, in this case, /i cannot be represented 
in a k-independent form as in the case of propagating waves. The physical explanation reads 
that since the EW does not have a Poynting vector in the direction of decay, it is the same for 
both the regular and the backward wave. Therefore, in this direction, there is no difference 
between LHM and RM. We have verified, however, (see Fig. 6) that if an EW propagates 
(has a real k) in a direction not parallel to the surface of the PC slab, it exhibits a negative 
refraction, similar to that shown in Fig. 4. 

Finally, we present a direct proof that the e and /i that describe the EW's are k-dependent 
and they are positive at the same frequency where e and \i for propagating modes are negative 
and k-independent. 
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First, we make two additional comments about Fig. 5(a): 

1) The plot Fig. 5(a) is completely the same if the slab of the RM is embedded into the 
HLHM and the materials are matched, i.e., they have the same absolute values of e and p. 

2) If the interfaces are not matched, the "reflected EW" would appear and there would be 
a maximum at the left boundary of the HLHM slab. 

To simplify the picture, we consider a PC slab with a BH surface so that surface waves are 
absent. The idea is to find e' and p! of the PC slab that are responsible for the propagation 
of EW's by choosing e and p of the homogenious regions in such a way that "reflected EW" 
is absent. 

Fig. 7 shows propagation of EW's through the PC slab surrounded by a HLHM (Fig. 
7(a)) and RM (Fig. 7(b)). To avoid a maximum at the left boundary of the PC slab in Fig. 
7(a) the HLHM should have e' = —1.78 and p! = —0.051. After that, one can see that the 
only difference between Figs. 7(a) and 5(a) is oscillations corresponding to the period of the 
PC. Thus, we get that e' = 1.78 and p! = 0.051 are responsible for the propagation of the 
EW's at a given k y and k through the PC. Fig. 7(b) shows propagation of a EW through 
the PC slab surrounded by a RM with e' = 1.78 and p! = 0.051. One can see that the EW 
decays inside the PC as exp(— kx) with the same value of k as in the RM. 

It follows from Fig. 7 that PC slab behaves as a RM with positive values of e' and p! at 
the same working frequency where propagating modes have negative e and p. Moreover, we 
found that if the values of k and k y of the EW are changed, keeping the same frequency, 
the values of e' and p! should be also changed. However, we still retain e'p' = n 2 . Thus, the 
parameters of the PC responsible for propagation of EW's are k-dependent. 



IV. IMAGE OF THE VESELAGO LENS 



Now we consider the image of the Veselago lens. The point source, as provided by our 
software, is the 2-D Green function H s = iH HQ 1 \pk ) exp — iutt, where k = uon/c and the 
Hankel function = J + iN . It is located in the RM at a distance a to the left of the 
PC slab. Similar to the geometry of Fig. 5, the distance a is 1/2 of the thickness of the PC 
slab. We simulated the Veselago lens with the same PC slab with AH cut. The parameters 
of the RM around the slab are obtained for the propagating modes and are given by Eq. 
(9). 
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FIG. 7: Magnetic field of the evanescent wave with k y = (\/5/2)feo, /c = 0.5&0j = un/c is shown 
with BH surface for the PC slab, a) The PC slab of the width lOd is surrounded by HLHM. b) The 
same PC slab is surrounded by regular media. The dashed line is a plot of the function exp (— kx). 



First, we show that the near-field image obtained with a thin lens is sharper than the far- 
field image of a thick lens (see Fig. 8). While the far- field image coincides with diffraction 
theory, the near-field image is beyond the diffraction limit. This is definitely connected 
with the amplification of EW's by the surface modes, as has been shown in the previous 
section and, following Pendry, this effect can be called "superlensing" . Fig. 8 shows the 
distributions of magnetic energy near the focus of both a thick lens and a thin lens in the 
lateral (y) direction. The distribution is symmetric for y and — y so it is shown for positive 
y only. Fig. 9 shows the distributions of magnetic energy near the focus of the thick lens in 
both lateral (y) and perpendicular (x) directions. 

Now we discuss the far-field image. The results of Sec. 3 give a possibility to calculate 
analytically the distribution of energy near the far-field focus of a thick Veselago lens with 
both AH and BH surfaces or of a thin lens with BH surface. In all these cases, evanescent 
waves do not reach the focus. At positive x, the Green function of the source can be 
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FIG. 8: Distribution of magnetic energy near the focus both for the thin lens and thick lens with 
AH surfaces in the lateral direction. The dashed line is the distribution of magnetic energy in the 
lateral direction for a thin lens; the solid line is the same for a thick lens. Analytical results as 
obtained from Eq. (13) for a thick lens are shown by circles. The dotted line shows the distribution 
of magnetic energy of the point source. The inset depicts the geometric meaning of k m . 

represented in a form H s = H p + H ev , where 

/fco exp i ( ky + x\/k\ — k 2 — out J 
7W=^ dk ( 10 ) 
-fc a/ k — k 

contains only propagating modes while 

r exp ( iky — x^ k 2 — k\ — iut J 

H ev = (H /n) / ^ == Uk (11) 

J |fc|>fco v — o 

contains only EW's. 

If the EW's decay inside the PC slab, only term H p contributes to the focus if the slab 
is thick enough. On the other hand, our results show that, for the propagating modes, the 
PC slab works exactly like a HLHM. Using the Fresnel equation, it is easy to show 3 that 
the field near the focus Hf is given by 

/fco exp i (ky + x'^/k 2 , — k 2 — ut ) 
jan? — Ldk - < 12 > 

Near the focus, the function Hf exactly coincides with Eq. (6) obtained in Sec. 2 by 
diffraction theory. Since diffraction theory in the far-field region is exact, it follows that 
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FIG. 9: Distribution of magnetic energy near the focus for the thick lens with AH surfaces in 
lateral and perpendicular directions. The solid line is the distribution of magnetic energy in the 
perpendicular direction; the dashed line is the same for the lateral direction. Circles show the 
analytical result for perpendicular direction, triangles show the same for lateral direction. Dotted 
line shows the distribution of magnetic energy of the point source. Analytical results are obtained 
fromEq. (13). 

EW's do not contribute to the far-field focus of Veselago lens independent of the microscopic 
nature of the LHM. 

All these results are obtained for the infinite lens. In our simulation, however, the slab 
has a finite length h in the y-direction. In this case, some propagating modes do not enter 
the slab. To take into account the finite aperture, one should integrate from — k m to k m in 
Eq. (12), where k m = k (h/2a)/^l + (h/2a) 2 and a is a distance from a source to a slab 
(see the inset of Fig. 8). Thus, 



In our case h/2a = 4,8 for thick and thin lenses, respectively. Due to this factor, the 
lateral distribution slightly differs from the square of the Bessel function Jq that follows 
from Eq. (7) for the infinite aperture. The results obtained from Eq. (13) are shown in 
Fig. 8 for thick (20d) lens. For the thick lens (far-field imaging), analytical calculations are 
in very good agreement with the computational data. Thus, Eq. (12) sets the diffraction 




k m exp % iky + x' 'a/^o — ^ 2 ~~ u ^ 



(13) 



dk. 
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FIG. 10: Simulation results for the distribution of magnetic energy near the focus in the lateral 
direction for two lenses with AH surfaces with different heights. The solid line is the distribution 
for the lens with h = 80d and the dashed line is for the lens with h = 120d. The dotted line is the 
plot of Jq (27ry / 'A) . The analytical results are almost indistinguishable from the simulation results. 
Llfigl0.eps 

limit for far-field imaging by a Veselago lens. In Fig. 10, we show the distribution along the 
lateral direction of two lenses with the same thickness a = 20d but a different h of 120c? and 
80c? correspondingly. As h becomes larger, the normalized intensity distribution approaches 
the square of the Bessel function. 

The first root of the Bessel function is at y = 0.38A. Thus, the sharp first peak of the 
Bessel-like function in the lateral direction can be confused with "superlensing" . We think, 
however, that any result that follows from Eq. (12) has nothing to do with "superlensing". 
Following Pendry, we would define superlensing as a result of amplification of EW's, but Eq. 
(10) does not contain EW's at all. It gives just the diffraction pattern of the Veselago lens 
that was first obtained in Ref. 5 for a 3-D lens. Then the energy distribution given by this 
equation should be called "regular lensing". 

The superposition principle is valid for combination of point sources with a different 
location. It can be used to create a focus with one maximum. To demonstrate this, we put 
three point sources along y axis with interval 0.3A successively, with the middle point source 
located at y — 0. In Fig. 11, we show both computational and analytical distribution of 
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FIG. 11: Distribution of magnetic energy along the lateral direction near the focus of three point 
sources with interval 0.3A with three point sources located at y = 0, y = — 0.3A and y = 0.3A, 
respectively (solid line) and the same for the interval 0.5A with three point sources located at 
y = 0, y = — 0.5A and y = 0.5A, respectively (dashed line). The analytical results are almost 
indistinguishable from the simulation results. 




FIG. 12: Distribution of magnetic energy of three point sources located along the y-axis with 
interval 0.3A. The distance between the sources and the PC slab is 1/2 of the slab thickness. 
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energy along the lateral direction near the focus and find that the image contains only one 
maximum value with the focal point at y — 0. The same picture shows that three point 
sources with separation 0.5A are already distinguishable by the Veselago lens. Fig. 12 shows 
the spatial distribution of the field intensity through all the system in the case of three point 
sources with interval 0.3A. 

An interesting feature of the far-field imaging follows from the above consideration. If 
the distribution of the object's field does not contain EW's, its far-field image created by 
an infinite Veselago lens will be perfect. The class of such field distributions includes all 
functions that can be represented in a form 

H(x,y) = J ° F(k) exp i (^ky + x^Jk% - k 2 - utj dk, (14) 

where F(k) is any integrable function in the interval [—ko,k }. The function H(x,y) obeys 
the homogeneous Helmholtz equation. Examples of such functions that produce perfect 
lateral imaging at x = are sm(k y)/y, \y\~ n J n (k \y\) and many others. 



V. CONCLUSION 



In conclusion, we show that a perfect image of the point source is impossible in both 
near-field and far-field regions and that the Veselago lens does not provide any superlensing 
in the far-field region. Both statements are independent of the nature of the LHM. We 
show that the far-field image has some peculiar features, such as perfect imaging for a whole 
class of functions. We present an analytical expression for the far-field imaging that sets 
a diffraction limit for the Veselago lens. We present an extensive study of the evanescent 
waves in a photonic crystal showing that a description of EW's can be done in the framework 
of a theory that includes the spatial dispersion of e and ji. As a result, both the universal 
amplification of EW's and the surface polaritons do not exist in the PC. However, in the 
near-field region (thin lenses) the amplification is possible if the surface supports regular 
surface waves. For the PC under study a thin lens with AH surface gives a substantially 
sharper image in the near- field region than that which follows from our diffraction limit. This 
is due to amplification of EW's shown in Fig. 5(c) and, following Pendry, this effect can be 
considered as a "superlensing" . However, the physics of this superlensing is connected with 
the properties of the surface rather than with the left-handed properties of the medium. 
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